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Abstract 

Ultra-high energy cosmic neutrinos are incisive probes of both astrophysical sources and new 
TeV-scale physics. Such neutrinos would create extensive air showers deep in the atmosphere. 
The absence of such showers implies upper limits on incoming neutrino fluxes and cross sections. 
Combining the exposures of AGASA, the largest existing ground array, with the exposure of the 
Fly's Eye fluorescence detector integrated over all its operating epochs, we derive 95% CL bounds 
that substantially improve existing limits. We begin with model-independent bounds on astro- 
physical fluxes, assuming standard model cross sections, and model- independent bounds on new 
physics cross sections, assuming a conservative cosmogenic flux. We then derive model-dependent 
constraints on new components of neutrino flux for several assumed power spectra, and we update 
bounds on the fundamental Planck scale in extra dimension scenarios from black hole produc- 
tion. For large numbers of extra dimensions, we find Mb > 2.0 (1.1) TeV for M§^ = Mb (5Mb), 
comparable to or exceeding the most stringent constraints to date. 
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I. INTRODUCTION 



Cosmic neutrinos provide a unique window on astrophysical processes because they escape 
from dense regions and typically propagate to the Earth unhindered 0. At ultra-high 
energies, they also provide an important probe of new ideas in particle physics. In contrast 
to all other standard model (SM) particles, their known interactions are so weak that new 
physics may easily alter neutrino properties, sometimes drastically. This is especially relevant 
for neutrinos with energies above 10 6 GeV, which interact with nucleons with center-of-mass 
energies above 1 TeV, where the SM is expected to be modified by new physics. 

The signal for ultra-high energy neutrinos is quasi-horizontal giant air showers initiated 
deep in the atmosphere [[J. This signal is well-studied and easily differentiated from air 
showers initiated by hadrons. The Earth's atmospheric depth rises from about 1000 g/cm 2 
vertically to nearly 36000 g/cm 2 horizontally. For all but the most extreme (and typically 
problematic || neutrino cross sections, the mean free path of neutrinos is larger than 
even the horizontal atmospheric depth. Neutrinos therefore interact with roughly equal 
probability at any point in the atmosphere and may initiate showers in the volume of air 
immediately above the detector. These will appear as "normal" showers, with large elec- 
tromagnetic components, curved fronts (a curvature radius of a few km), and signals well 
spread over time (of the order of microseconds). 

On the other hand, hadrons have interaction lengths of the order of 40 g/cm 2 and so 
always interact at the top of the atmosphere. The electromagnetic component of an air 
shower has mean interaction length ~ 45 — 60 g/cm 2 . For a quasi-horizontal shower initiated 
by an ordinary hadron, then, this component is absorbed long before reaching the ground, 
as it has passed through the equivalent of several vertical atmospheres — 2 at a zenith angle 
of 60°, 3 at 70°, and 6 at 80°. In these showers, only high energy muons created in the first 
few generations of particles survive past 2 equivalent vertical atmospheres. The shape of 
the resulting shower front is therefore very flat (with curvature radius above 100 km), and 
its time extension is very short (less than 50 ns). 

These shower characteristics are exploited by both ground arrays and fluorescence detec- 
tors in searches for primary cosmic ray neutrinos. At present, no ultra-high energy neutrino 
signal has been reported. Here we determine the total exposure for neutrino detection 
from existing facilities and derive both model-independent and model-dependent bounds on 
astrophysical neutrino fluxes and new neutrino interactions. 

The outline of the paper is as follows. In Sec. [TT] we examine acceptances for neutrino 
detection and compute the current combined total exposure using all available data from 
the Akeno Giant Air Shower Array (AGASA) || and Fly's Eye || experiments. In Sec. Il l 
we determine model-independent bounds on the total neutrino flux, assuming SM cross 
sections. To derive model-independent results, we assume only that fluxes are confined to a 
small window around some central neutrino energy and obtain bounds as a function of this 
central energy. After that, we assume a power law neutrino flux d<&/dE v oc E^ 1 to obtain 
stronger, but more model-dependent, bounds on the total neutrino flux from integrating 
over all energies. 

In Sec. |y| we derive model-independent bounds on high energy neutrino cross sections, 
assuming a conservative cosmogenic flux. These significantly improve existing limits Jj]]. We 
then derive model-dependent bounds on cross sections, focusing on the example of TeV-scale 
gravity scenarios in Sec. |V|. We improve existing constraints [H on the fundamental Planck 
scale from the non-observation of microscopic black hole production by cosmic neutrinos 0, 
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10, [11], [12|, [H|. For large numbers of extra dimensions, these bounds are comparable to or 
exceed all existing bounds on extra dimensions. Sec. |VI| contains our conclusions. 



II. NEUTRINO EXPOSURE 



To estimate the sensitivity of a ground-based detector to neutrino-initiated showers, we 
must first compute its exposure for various types of showers. The exposure is the product 
of the effective aperture and the range of depths within which the shower must originate to 
trigger the device, integrated over time. The effective aperture is the detector's projected 
area weighted by detection probability and integrated over solid angle. The exposure is, 
then, 

£(E sh ) « f T dt (AQ) cS (E sh , t) e~ h l H dh , (1) 

JO JO Pwater 

where T is the total observation time of the detector, /i max = 15 km, H « 8 km, and 

Po ~ 1.15 x 10~ 3 p wat er is the density of the atmosphere at ground level. The effective 
aperture is Jl4]] 

(An) cS (E sh , t) = / A{t) V(E sh , 6, t) 2vr sin 6 d9 (2) 

$min 

where A(t) is the detector's area, V(E B i x ,9,t) is the probability that a shower that arrives 
with energy E s ^ and zenith angle 8 triggers the detector, and the angular cuts 9 min and # max 
are chosen to optimize detection efficiency while eliminating hadronic background. 

Exposures at a given detector depend on detection method and shower type. We will 
refer to showers as hadronic or electromagnetic (EM), depending on the nature of their 
first interaction, irrespective of their later development. With this convention, for example, 
the SM neutral current process vN — > vX produces a hadronic shower, while the charged 
current process for electron neutrinos i/ e N — > eX produces both a hadronic shower and an 
EM shower. 



For ground arrays, exposures for hadronic and EM showers differ [15|. In hadronic show- 



ers, the initial hadronic interaction produces a strong muon component that remains until 
the shower reaches the ground. This muon component is largely absent for EM showers. 
These muons significantly enhance triggering efficiencies for ground arrays, which are sensi- 
tive only to ground level activity, and so exposures for hadronic showers exceed exposures 
for EM showers at ground arrays. Above some critical energy, which depends on the total 
effective area of the array, the detector exposure saturates. Exposures for hadronic and 
EM showers for the AGASA ground array || are given in Fig. [I]. These exposures are for 
1.5 x 10 8 s of livetime between December 1995 and November 2000. The hadronic exposure 
was derived in Ref. ||, based on results from searches for deeply penetrating showers 



and conservative assumptions. We derive the EM exposure by comparison with results for 
the Auger experiment. Effective apertures for Auger have been calculated in Ref. [15]. For 
AGASA, we adopt the Auger aperture for quasi-horizontal EM showers with axes falling 
in the array, reduced by a factor of 30, the ratio of surface areas of the arrays. We have 



checked that, using this EM exposure, we reproduce AGASA's bounds |17j on v e fluxes to 
within 20%. 

In contrast to ground arrays, the hadronic and EM exposures of fluorescence detectors 
are very similar. Fluorescence detectors are sensitive to the total EM activity along the 
entire longitudinal development of the shower. In EM showers, essentially all of the energy 
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FIG. 1: The figure shows the total exposures for hadronic showers at AGASA (solid), EM showers 
at AGASA (dashed), and EM and hadronic showers at Fly's Eye (dotted) as functions of shower 



energy. 



produces EM activity. In hadronic showers, hadronic collisions produce equal numbers of 7r°, 
7r + , and tt~ . The 7r° decays to photons, producing EM energy. However, the charged pions 
typically interact before decaying. Through successive interactions, most of their energy 
also becomes EM activity. As a result, roughly 90% of the energy in hadronic showers is 
EM. The response and efficiency of fluorescence detectors to hadronic and EM showers is 
therefore expected to be similar ]IB |. In this work we adopt the total Fly's Eye exposure 



reported in Ref. |J for both hadronic and EM showers. This exposure, given in Fig. [I], 
includes not only data from the first epoch of observation (February 1983 to May 1985), 
most of which were reported in Ref. [|18| , but also data from four additional running periods 
between November 1985 to July 1992. The additional periods enhance the total Fly's Eye 
exposure by roughly a factor of 3. 

The AGASA Collaboration has searched for quasi-horizontal showers that are deeply- 
penetrating, with depth at shower maximum X max > 2500 g/cm 2 [17|. At AGASA, the 
location of the shower maximum is determined through its correlation to two measurable 
quantities: rj, which parameterizes the lateral distribution of charged particles at ground 
level, and 5, which parameterizes the curvature of the shower front. Deeply penetrating 
events must satisfy X^ ax , X^ ax > 2500 g/cm 2 . The expected background from hadronic 
showers is I^Iq^Iq;^, where the first uncertainty is from Monte Carlo statistics, and the 
second is systematic. Among the 6 candidate events, 5 have values of and/or X^ ax 

that barely exceed 2500 g/cm 2 , and are well within AX max of this value, where AX max is 
the estimated precision with which X max can be reconstructed. The AGASA Collaboration 
thus concludes that there is no significant enhancement of deeply penetrating shower rates 
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given the detector's resolution. 

The Fly's Eye observes an air shower as a nitrogen fluorescence light source traveling 
at the speed of light through the atmosphere The light is emitted isotropically with 
intensity proportional to the number of charged particles in the shower. The received light 
profile is reconstructed by a 3 parameter fit to the charged particle density (largely electrons 
and positrons) along the shower track. The parameters are the depth of the observed first 
interaction X ; the depth of the shower maximum X max ; and the density A max at X max . In 
a running time of eleven years, Fly's Eye recorded more than 5000 events. However, there 
are no neutrino candidates with shower maximum deeper than 2500 g/cm 2 [T8], |T5|| . 

All in all, given 1 event that unambiguously passes all cuts with 1.72 events expected 
from hadronic background, the combined results of AGASA and Fly's Eye imply an upper 
bound of 3.5 events at 95% CL IpOfl from neutrino fluxes. 



The event rate for quasi-horizontal deep showers from ultra-high energy neutrinos is 

N = Y, fdE.NA^- a^xiEi) , (3) 

i,X J dh i 

where the sum is over all neutrino species % = v e , v e , z/ M , z/ M , u T , v T , and all final states X. 
Na = 6.022 x 10 23 is Avogadro's number, and d^i/dEi is the source flux of neutrino species 
i. Finally, £ix{Ei) is the appropriate exposure measured in cm 3 we sr • time. 

To clarify, we present appropriate exposures for some SM processes. At the ultra-high 
energies of interest, on average 20% of the total neutrino energy goes into hadronic recoil for 
both SM neutral current and charged current events ||21|| . Exposures for SM charged current 
events at AGASA are therefore 

Su e x(Eu e ) = min{£ had (0.2£j+£ EM (0.8£j,£ sat } (4) 
£uM E »J = 4ad(0.2^J (5) 
£u t x{E Vt ) = £ had (0.2E„ T ) , (6) 

with identical expressions for anti-neutrinos. The exposures for and v T are identical, 
because at these energies, r leptons typically do not decay before arriving at the Earth's 
surface. For z/ e , the expression includes the effects of saturation, noted previously. For 
AGASA, the exposure saturates at £ sat ~ 5.3 km 3 we sr yr at E s ^ m 10 10 GeV (see Fig. [j]). 
For neutral current interactions at AGASA, all of the right-hand entries are £^had(0.2£' ;/i ). 
For Fly's Eye, the corresponding charged current exposures are simply 

£ Ve x{E Ve ) = £{E Ue ) (7) 

£uM E ^) = £(0-2£„J (8) 

£ VtX {E Vt ) = £{Q.2E Vt ) . (9) 

For neutral current exposures at Fly's Eye, all of the right-hand entries are £(0.2i^.). 

Given the exposures £ described above, the absence of events implies upper bounds on 
cross sections, assuming some fixed flux, or upper bounds on fluxes, assuming some fixed 
cross section. We now consider these two possibilities in turn. 



III. BOUNDS ON ASTROPHYSICAL NEUTRINO FLUXES 

To derive bounds on the neutrino flux, we assume SM charged and neutral current inter- 
actions for all neutrinos. We further assume that the source flux of neutrinos, which in the 



5 



energy range of interest is dominantly and u e at production, is completely mixed in 

flavor upon arrival at the Earth. There is now strong evidence for maximal mixing among all 
neutrino species ||22|| . In addition, given mass differences of Am 2 > 10~ 6 eV 2 , the neutrino 
flux is expected to be completely mixed if they travel a distance > 0.1 Mpc. The assumption 
of equal flavor representation upon arrival is therefore strongly supported by data |£J . 

We first derive model-independent bounds on the total neutrino flux. Let us start by 
noting that if the number of events integrated over energy is bounded by 3.5, then it is 
certainly true bin by bin in energy. Thus, using Eq. ([5J) we obtain 

]T / dEiNA-^aiN^xiE^SixiEi) < 3.5 , (10) 

at 95% CL for some interval A. Here the sum over X takes into account charge and neutral 
current processes. In a logarithmic interval A where a single power law approximation 

d® 

=-± (TiN-xiEk) EuciEi) ~ £f (11) 
at,i 

is valid, a straightforward calculation shows that 

r {E)e A / 2 (IE- d$- smb. 5 

/ / -=± Ei — f (T iN ^ x {Ei) Six(Ei) = (a ViN ^ x (Ei) £ iX (Ei) E { d^/dE^^—A , 

J(E)e- A / 2 Ei dEi 5 

(12) 

where 5 = (a + l)A/2 and (A) denotes the quantity A evaluated at the center of the 
logarithmic interval. The parameter a = 0.363 + /3 — 7, where the 0.363 is the power law 
index of the SM neutrino cross sections ]2T] , and f3 and —7 are the power law indices (in the 
interval A) of the exposure and flux d^i/dEi, respectively. Since sinh 5/5 > 1, a conservative 
bound may be obtained from Eqs. ([IT]) and (|12|) : 

N A Y,(°nN->x{Ei)) (£ix(Ei)) (EidQi/dEi) < 3.5/A . (13) 

i,X 

In this work we choose A = 1, corresponding to one e-folding of energy, as a likely interval 
in which the single power law behavior is valid. By setting (E^d^i/dEi) = ^(E u d& u /dE u ), 
where is the total neutrino flux, we obtain model-independent upper limits on the total 
neutrino flux at 95% CL. The results are given in Table |. 

These model-independent upper bounds on the total neutrino flux can be strengthened 
by assuming a particular flux behavior. For example, if the neutrino flux falls like 



dE„ V Ei 



Eq. ([TOD leads to J < 1.4 x KT 5 km" 2 sr" 1 yr" 1 GeV" 1 , for 7 = 2 and E = 10 s GeV. Under 
the same assumptions for 7 = 1.5 one obtains Jo < 9.8 x 10 -7 km -2 sr" 1 yr" 1 GeV" 1 . Fig. |2| 
shows both the model-independent bounds of Table | as well as the bounds on flux under 
the power law assumptions just discussed. Additionally displayed in the figure are the upper 
limits on the + v e flux obtained from the non-observation of microwave Cerenkov pulses 
from EM showers induced by neutrinos in the Moon's rim, as measured by the Goldstone 
Lunar Ultra-high energy neutrino Experiment (GLUE) [p4[| , as well as bounds from the 
search for radio pulses from EM showers created by electron neutrino collisions in ice at 
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TABLE I: Model-independent upper limits on the differential neutrino flux at 95% CL. 



E v (GeV) {E u d<5> v /dE u ) (km- 2 sr" 1 yr" 1 ) 

1 x 10 8 1.8 x 10 5 

3 x 10 8 4.1 x 10 4 

1 x 10 9 7.9 x 10 3 

3 x 10 9 2.2 x 10 3 

1 x 10 10 5.0 x 10 2 

3 x 10 10 1.6 x 10 2 

1 x 10 11 6.8 x 10 1 



the Radio Ice Cerenkov Experiment (RICE) |[25|| . Comparing model- independent bounds to 
model-independent bounds, and model-dependent bounds to model-dependent bounds, we 
find that the bounds obtained here are significantly more stringent than existing limits. 

Also displayed in Fig. [2] as a single point with error bars is the total neutrino flux required 
by Z-bursts f28j. The Z-burst model proposes that the observed ultra-high energy cosmic 
rays with energies above 10 11 GeV are secondaries resulting from resonant annihilation of 
ultra-high energy neutrinos on relic neutrinos at the Z pole ||29fl . A recent analysis f27 



includes several possibilities for a diffuse background of protons, as distinct from protons 
resulting from the Z-burst itself. The point shown in the figure corresponds to the case in 
which the proton background originates at distances < 50 Mpc. The normalized Hubble 
expansion rate, the matter and vacuum energy densities, and the maximum redshift are 
taken as h — 0.71, Qm = 0.3, Qa = 0.7, and z max = 2, respectively. The neutrino flux 
is assumed to have no cosmological evolution. The horizontal errors result from the la 
uncertainty in the neutrino mass determination. The errors in the flux reflect the statistical 
fluctuations in the fits, as well as the uncertainty in the Hubble expansion rate. 

A more speculative explanation of the mysterious events at the high end of the spectrum 
assumes that the cosmic ray primaries arise in the decay of massive elementary X particles. 
Sources of these exotic particles could be either topological defects left over from early 
universe phase transitions (mx ~ 10 16 — 10 19 GeV) ||30|| , or some long-lived metastable 
superheavy {m x ^ 10 12 GeV) relic particles produced through vacuum fluctuations during 
the inflationary stage of the universe pfl] . The X particles typically decay to leptons and 
quarks. The latter produce jets of hadrons containing mainly pions, together with a 3% 
admixture of nucleons. The predicted spectrum would thus be dominated by gamma rays 
and neutrinos produced via pion decay. The neutrino flux bounds derived in this work 
therefore seriously constrain this type of model. Moreover, recent analyses of Haverah 
Park data [|32| suggest that less than 50% of the primary cosmic rays above 4 x 10 10 GeV 
can be photons at 95% CL. Note that mechanisms which successfully deplete the high 
energy photons (such as efficient absorption on the universal and galactic radio background) 
require an increase in the neutrino flux to maintain the overall normalization of the observed 
spectrum | 33fl . Definite quantitative comparison with the neutrino flux bounds presented 
here can be obtained by specifying the nature of the decay of the X particle. 
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FIG. 2: 95% CL model-independent (squares) and model-dependent (solid lines) upper limits on 
the total neutrino flux E u dQ u /dE v derived here using the combined exposures of AGASA and Fly's 
Eye. The model-dependent bounds assume d§ u /dE u oc E~^, with 7 = 1.5 and 2. For comparison, 
we also present limits on the + v e from GLUE (circles) [^4[ and the 95% CL bound from RICE 
on the v e flux assuming a power law spectrum with 7 = 2 (dashed) p5| . The dash-dotted line is 
the predicted total cosmogenic flux from pion photo-production |^], and the point with error bars 
is the total neutrino flux required by Z-burst models, as derived in p7|]. 



IV. BOUNDS ON NEW PHYSICS INTERACTIONS 



In this section we examine the potential of probing the big desert that lies between the 
electroweak and grand unified theory (GUT) scales using ultra-high energy neutrino interac- 
tions. To derive bounds on possible new physics contributions to neutrino cross sections, we 
assume the "guaranteed" flux of cosmogenic neutrinos arising from pion photo-production 
from ultra-high energy protons propagating through the cosmic microwave background. This 
flux depends on the cosmological evolution of the cosmic ray sources. Throughout this work, 
we conservatively adopt the estimates of Protheroe and Johnson [26] with nucleon source 
spectrum scaling as d§ nuclcon /dE oc E~ 2 and extending up to the cutoff energy 10 12,5 GeV. 
We assume also a cosmological source evolution scaling as (1 + z) A for redshift z < 1.9 fl34 . 
The total ultra-high energy cosmogenic neutrino flux is shown in Fig. |2|. 

We consider this flux highly conservative. For example, one might conjecture that the 
cosmogenic flux is absent because the observed cosmic rays with energies > 10 11 GeV are 
protons or nuclei generated by nearby sources within 50 Mpc. However, none of the known 
nearby candidate sources, such as Virgo, M82, Centaurus A, and galactic pulsars [^] is as 
powerful as more distant sources, such as Cygnus A and Pictor A ||36|| . The latter must 
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TABLE II: 95% CL model-independent upper limits on the neutrino-nucleon cross section. 



E u (GeV) 


°In (P b ) 


°In A (P b ) 


1 x 10 10 


1.8 x 10 6 


7.4 x 10 6 


3 x 10 10 


8.8 x 10 6 


2.2 x 10 7 


1 x 10 11 


8.1 x 10 7 


1.1 x 10 8 



therefore inject nucleons with energies of at least 10 12 GeV. Indeed, contributions to the 
nucleon channel from semi-local sources p7| , as well as recently discussed possibilities |38 



of source spectra harder than E~ 1,5 and source evolutions stronger than (1 + z) 4 , would all 
significantly enhance the cosmogenic neutrino flux. Contributions from decaying topological 
defects, active galactic nuclei, and other speculative sources would have a similar effect. If 
realized in nature, any one of these possibilities would strengthen our bounds. 

For SM interactions and the cosmogenic neutrino flux given in Fig. ^|, the expected rate 
for deeply penetrating showers at AGASA and Fly's Eye is about 0.02 events per year, 
and so negligible. New physics may contribute to neutrino cross sections in a multitude of 
ways || [11], [39], 40], with different contributions for different flavors, and a variety of final 



states producing showers with a variety of hadronic and EM shower components. Here we 
assume that the new physics is flavor-blind, inducing equivalent cross section modifications 
for all neutrino species, and that the resulting new physics final state leads to showers with 
negligible EM component. We then consider two simple but representative cases: (1) y = 1, 
and (2) y = 0.1, where y = E sh / E v is the average inelasticity. In case (1), the shower energy 
-Esh differs little from the neutrino energy E U) as in TeV-scale black hole production (see 
below). In case (2), E^ is substantially less than E u . This holds, for example, when neutral 
current interactions are enhanced by the exchange of Kaluza-Klein (KK) gravitons || . 

Our starting point to derive model-independent bounds on the total neutrino-nucleon 
cross section a U N is again Eq. (|13|). For reasons given previously, we choose A = 1, so that 
Eq. ( |13|) becomes 



N A (a vN ^ x (E v )) (S(yE v )) (E v d<$> v / dE v ) < 3.5 . (15) 

Using the exposures in Fig. [I] and the cosmogenic flux in Fig. [|, we find mo del- independent 
bounds on the neutrino cross section for different energies and the two inelasticities above. 
The resulting limits are given in Table |I] and shown in Fig. |3|. For reference, the SM cross 



sections [^TJ are also given in the figure. 

These bounds assume that neutrinos produce deeply penetrating quasi-horizontal show- 
ers. They may be avoided if neutrinos are so strongly interacting that they shower high in the 
atmosphere. For zenith angles of 70°, this requires an interaction length below 3000 g/cm 2 , 
corresponding to a cross section above 6 x 10 8 pb. If the interaction length is between 3000 
g/cm 2 and the horizontal atmospheric depth of 36000 g/cm 2 , corresponding to a vN cross 
section between 6 x 10 8 pb and 5 x 10 7 pb, respectively, an exponential decrease in the 
event rate for showers of a given total energy will be observed beyond a critical angle. Since 
this range includes our upper bounds on a U N at 10 11 GeV, it is conceivable that <j v n could 
thus be measured directly. Bounds on anomalous vN cross sections which are near-hadronic 
strength can be obtained through considerations of their feedback to low energy physics via 



dispersion relations [[41 
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FIG. 3: 95% CL model-independent upper limits on a v N for inelasticity y = 1 (filled squares) and 
y = 0.1 (open squares). The solid contour is the upper limit on the black hole production cross 
section for n = 7 and x m i n = 5, corresponding to Mjj = 1.1 TeV. For comparison, the SM neutral 
current, charged current, and total neutrino nucleon cross sections are also given (dot-dashed, from 
below). These bounds assume that the neutrino cross sections does not exceed 6 x 10 8 pb (see 
text). 



The bounds shown in Table [TT] strengthen previous bounds |7J by roughly one order of 
magnitude. This enhancement results from a number of factors. First, the exposure used 
in Ref. |7j was limited to the first 6 x 10 6 s of Fly's Eye operation. The updated exposure 
used here is roughly 9 times larger. Second, the cosmogenic flux used in Ref. |7j assumes a 
source cutoff energy of 10 11 ' 5 GeV, and is smaller by roughly a factor of 4 than the one used 
here. We find it unnatural to assume that the the source cutoff coincides with the maximum 
observed cosmic ray energy, and so have used the flux corresponding to the larger cutoff 
energy. Finally, we have presented 95% CL limits, corresponding to a limit of 3.5 events. 

As in the case of flux bounds, the model-independent cross section bounds of Table |Tj 
will be improved if one assumes a cross section shape and so can integrate over all energies. 
We consider a particular example in the next section. 



V. IMPLICATIONS FOR TEV-SCALE GRAVITY 

The idea that our universe could be a brane embedded in some higher dimensional world 



has received a great deal of renewed attention over the last 5 years E2, 43 . From a phe- 



nomenological point of view, this possibility presents a new perspective on the hierarchy 
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between the gravitational and electroweak mass scales. In these scenarios, the effective 
4- dimensional Planck scale Mpi ~ 10 19 GeV is determined by the fundamental (4 + n)- 
dimensional Planck scale Mp ~ 1 TeV and the geometry of the n extra dimensions. 

Arguably the most fascinating prediction of TeV-scale gravity is the production of black 
holes (BHs) in observable collisions of elementary particles [13, For cosmic rays, this 
implies that ultra-high energy neutrinos may produce BHs in the atmosphere, initiating deep 
quasi- horizontal showers far above SM rate || . BH production therefore provides a specific 
example of a model-dependent cross section that is bounded by the arguments discussed 
above. 

TeV-scale gravity also has a number of other implications for cosmic rays. The implica- 
tions of perturbative KK graviton exchange has been considered in a number of scenarios J39J . 
However, for cosmic rays, in contrast to the case at colliders, there is an abundance of center- 
of-mass energy. Extra dimensional effects will therefore first appear as non-perturbative BH 
production in processes with center-of-mass energies above the fundamental Planck scale, 
rather than through perturbative effects below the Planck scale. This is in stark contrast to 
the case at colliders, where the sensitivity to KK graviton effects surpasses the sensitivity 
to BH production. 1 

The sensitivity of current cosmic ray experiments to BH production, as well as that of 
facilities expected in the not-too-distant future, has been thoroughly investigated || |9], [11], 
12, 13| . The parton-parton cross section is estimated from the geometric area of the BH 

2 n 



horizon and is of order <j,- 



where 



,(Mj 



BH I 



D 



Mi 



1311 



M, 



D 



i 

l + n 



2™7r(™-3)/2r(n±3) 
n + 2 



i 

l + n 



(16) 



is the radius of a Schwarzschild BH in 4+n dimensions HGj. Criticisms of the absorptive 
black disc scattering amplitude, which center on the exponential suppression of transitions 
involving a (few-particle) quantum state to a (many-particle) semiclassical state j|7|], have 
been addressed in Refs. f|8, ^9|. The geometric cross section applies for both flat and 
hyperbolic |50] extra dimensions that are larger than the Schwarzschild radius, and for 
warped extra dimensions where r s is small compared to the curvature scale of the geometry 
associated with the warped subspace [[51], [52|| . 

The total production cross section for BHs with mass Mbh = 



sx is then 



°Viv— >bh ( E u ) = V / dx &i(Vxs) fi(x, Q) 



(17) 



where s = 2m^Ey^ x is the parton momentum fraction, the /j are parton distribution 
functions (pdf's), M^ 11 is the minimum BH mass, and the sum is carried out over all partons 
in the nucleon. The choice of the momentum transfer Q is governed by considering the time 
or distance scale probed by the interaction. According to Thome's hoop conjecture f53fl , the 
formation of a well-defined horizon in four dimensions occurs when the colliding particles 
are at a distance ~ r s apart. (Note that there could be an n-dependent factor for higher 
number of dimensions.) This has led to the advocacy of the choice Q ~ r~ l | 1T| , which 
has the advantage of a sensible limit at very high energies. However, as has been pointed 



We thank S. Dimopoulos for emphasizing this point. 
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out by Dimopoulos and Emparan p8 |, string progenitors can give experimental signals akin 



to BH for the collision energies and values of under present consideration. In cosmic 

ray experiments, these signals are indistinguishable from those of BH decay (more on this 
below). Detailed calculations |54| show that these "string ball" cross sections can exceed BH 
cross sections. In this region, the dual resonance picture of string theory would suggest a 
choice Q ~ M res ~ ^/xs. Since we have chosen to use only BH cross sections over the entire 
energy region, we set Q = min{MBH, 10 TeV}, where the upper limit is from the CTEQ5M1 
distribution functions We are aware that for Q > string scale, the pdf's will receive 
significant corrections from the rapid increase of degrees of freedom. Fortunately, as noted 
in Ref. |§, the cross section <j u n~>bk is largely insensitive to the details of the choice of Q. 
For example, the two choices discussed here result in cross sections that differ by only 10% 
to 20%. 

Once produced, BHs will Hawking evaporate with a temperature proportional to the 
inverse radius Tjj = (n + l)/(47rr s ). The wavelength A = 2tt/Th corresponding to this 
temperature is larger than the BH size. Hence, to first approximation the BH behaves like 
a point-radiator with entropy 

n + 2 K J 

and mean lifetime 

3+n 

T » ~ W D llfo J ' (19) 

Microscopic black holes therefore decay almost instantaneously to a thermal distribution of 
SM particles. As very few SM particles are invisible to cosmic ray detectors, the neutrino en- 
ergy is almost entirely transformed into shower energy. The EM component of these showers 
differs substantially from that of SM neutrino interactions, allowing a good characterization 
of the phenomenon against background when the BH entropy 3> 1 jllj. Recently, it has 



been noted that BH recoil induced by KK-graviton emission may launch the BH out of the 



brane ||56|| . In this case, BH radiation would be prematurely terminated from the perspec- 
tive of brane observers. This effect would drastically deplete the rate of deeply developing 
showers, and could be misinterpreted as a sharp cutoff on the ultra-high energy neutrino 
spectrum. We assume here that the effect of recoil is negligible. 

An important parameter in determining the BH cross section is x m ; n = M^/Mp, the 
ratio of the minimal black hole mass to the fundamental Planck scale. The above description 
of BH production and decay relies on semi-classical arguments, valid for large x m i n or, 
equivalently, large BH entropy. For large x m i n , thermal fluctuations due to particle emission 
are small (S ^> 1) [|57 |, statistical fluctuations in the microcanonical ensemble are small 



(VS 3> 1), and quantum gravity effects may be safely neglected. In addition, gravitational 
effects of the brane on BH production, which are ignored in all analyses to date, are expected 
to be insignificant for BH masses well above the brane tension, which is presumably ~ Mp. 
None of these is true for x min ~ 1. 

Cosmic ray experiments are largely insensitive to the exact details of BH decay. Whatever 
happens ~ 1, it seems quite reasonable to expect that BHs or their Planck mass 

progenitors will decay visibly, triggering deeply atmospheric cascade developments. There 
is, however, sensitivity to the BH production cross section. In string theory, BH production 
is expected to gradually pass to the regime of string ball production as M B h approaches M D . 
Evidence from this picture suggests that the BH production cross section is not radically 
altered in this limit [13]. This addresses many of the concerns listed above, but does not 
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FIG. 4: 95% CL lower limits on the fundamental Planck scale as a function of x m i n for n = 1, . . . , 7 
extra dimensions (from below). 



address the problem of brane effects on BH production 
« 1. In our analysis, we avoid choosing a specific x m 



x - 



these may still be large for 
rather, we present results for 
the generous range 1 < x m \ n < 10. As we will see, the bounds are rather insensitive to x m i n , 
in contrast to the case at colliders such as the LHC. 

In Fig. H we show the lower limit on Mr> as a function of x m ; n corresponding to < 3.5 
events (95% CL) observed for the combined exposures of AGASA and Fly's Eye. Bounds 
on Mjj from table-top gravity experiments, as well as from astrophysical and cosmological 
considerations, greatly exceed 1 TeV, for models with n < 4 flat extra dimensions [j5"5| . For 
n > 4, however, the bounds of Fig. ^ are among the most stringent to date. Note that these 
bounds do not depend on the shape of the extra dimensions and are valid for both warped 
and non-warped scenarios 



For x m i n = 1, the bounds extend up to 2.0 TeV for n = 7. 
Moreover, assuming x min = 3, for which the entropy S > 10, the bounds derived with the 
combined exposure, for n = 5,6,7, are M D > 1.26 TeV, 1.30 TeV, 1.40 TeV, respectively. 
All of these exceed bounds from the Tevatron and LEP ||59|| , even in the case where the 
brane softening parameter A |Kj is as large as M D (see || for details). 2 

In Fig. m we have also plotted the maximal BH cross section, corresponding to Mo = 
1.1 TeV, for the case n = 7 and x m i n = 5. As expected, given a model for the cross section's 
energy dependence, the resulting bounds on new interactions are much more stringent than 
the model-independent limits derived above. 



2 The increased exposure used here will also strengthen existing bounds |pl| from p-brane production |62[ 
in asymmetric compactifications. 
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FIG. 5: Black hole events at Auger in 3 years (1, 10, 100, 1000, from above) (solid) and at the 
LHC for integrated luminosity 100 fb^ 1 (1, 10, 100, ... , 10 , from above) (dashed) for n = 7 
extra dimensions. 



Finally, as noted above, the event rates for black hole production by cosmic rays are 
fairly insensitive to the choice of x m ; n . This contrasts sharply with the case at colliders. 
Specifically, the total production cross section of a BH of mass Mrh ~ \/Fs in a pp collision 
is given by 

Opp^BHOmin, s) = / d,T Mx)fj(r/x)&ij, (20) 



where r is the parton-parton center-of-mass energy squared fraction, and y/T min s is the 
minimum center-of-mass energy for which the black disc approximation is valid. The number 
of BH produced at the LHC (y/s = 14 TeV and luminosity C = 10 34 cm~ 2 s _1 fB5| ) is then 



N Lnc = J (Tpp^BH £> dt. In Fig. ^|, we show both Auger [Q and LHC event rates for various 
x m i n . For fixed Md, the LHC event rates drop by one to two orders of magnitude for 
every unit increase in x m i n , while the Auguer event rates are relatively stable. This may 
be understood as resulting from a combination of the very high energies available in cosmic 
neutrinos and the fact that the parton energy is not degraded by a parton distribution 
function in the cosmic neutrino 'beam.' 



VI. CONCLUSIONS 

In the first part of this paper, we derived new limits on the cosmic neutrino flux striking 
the Earth's atmosphere. This was accomplished by searching for quasi-horizontal deeply 
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developing showers in ultra-high energy cosmic ray data, taking into account the combined 
exposures of the AGASA and Fly's Eye experiments. Our results significantly strengthen 
existing limits and present serious problems for models where exotic elementary X particles 
cascade decay to cosmic ray particles. In particular, models where topological defects are 
responsible for the events detected with energies > 10 11 GeV are severely constrained, be- 
cause neutrinos are typically a significant component in X decays, and have a hard spectrum 
extending up to Mqut ~ 10 16 GeV. The bounds obtained in this paper will also challenge 
any attempt to normalize the observed spectrum to the proton flux as predicted by top 
down models. 

In the second part of the paper, we used the atmosphere as a giant calorimeter to probe 
neutrino-nucleon cross sections at -\/s > 1 TeV. We first combined the complete neutrino 
exposure of the above-mentioned facilities with the flux of cosmogenic neutrinos, to de- 
rive model-independent upper bounds on the neutrino-nucleon cross section. These bounds 
strengthen existing limits by roughly one order of magnitude. We then considered TeV- 
scale gravity models to study BH production. The upper bounds on the neutrino-nucleon 
cross section implied lower limits on the fundamental Planck scale, which represent the best 
existing limits on TeV-scale gravity for n > 5 extra spatial dimensions. 
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